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ABSTRACT 


The development and structure of Dean vortices ina 
curved channel were measured and studied. A 40:1 aspect 
ratio channel with mild curvature was employed. Circulation 
of mean streamwise vorticity, and maxima of streamwise, 
Spanwise and radial components of vorticity increase as the 
Dean number increases from 40 to 140. Mean vorticity 
components were determined from measurements of the three 
components of mean velocity using a miniature five hole 
pressure probe. Spectra obtained from a subminiature 
crossed-wire probe show principal and harmonic frequencies 
for Dean numbers greater than 165, when vorticies are 
twisting. Spectral peaks at other Dean numbers provide 
evidence of oscillations also observed during flow 
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I. INTRODUCTION 


A. BACKGROUND 

Flow in a curved channel with mild curvature and a 40 
to 1 aspect ratio was studied. Results are presented for 
channel Reynolds numbers, Ud/v , less than 1600. [In 
particular, the structure of Dean vortices and their 
development as Reynolds numbers increases was 
investigated. The presence of Dean vortices in curved 
channel flow has an important influence on flow structure, 
heat transfer characteristics and eventually, transition 
to turbulent flow. 

Centrifugal forces imposed by concave curvature of 
streamlines cause curved channel flow to be non-azimuthal. 
As curvature is imposed, centrifugal forces are greater on 
the faster moving fluid particles near the center of the 
channel than near the concave or convex surfaces. The 
centrifugal instability causes the faster moving fluid to 
displace fluid near the concave wall. The particles near 
the concave surface are then forced in the spanwise 
direction and then toward the center of the channel. This 
centrifugal instability results in spanwise periodic pairs 


of counter-rotating vortices (Figures 1, 2) [{Ref. 1]. 


While significant experimental and analytic research has 
been accomplished on flows between concentric cylinders and 
over concave surfaces, little work has been done to study 
flow in curved channels. The first analytical study of 
curved channel flow was by W.R. Dean [Ref. 2] in 1928. He 
found that, for a narrow gapped channel (r;/r, ~ 1), flow is 
unstable to small disturbances when the value of Ud/v /d/r; 
1s greater than 35.92. The parameter Rel ey d/o eemod 
known as the Dean number, De. Reid (Ref. 3) and Hammerlin 
(Ref. 4] verified Dean’s findings in additional analytic 
work and further defined a neutral stability curve which 
divides stable and unstable flow regimes. 

The first published experimental work done on high 
aspect ratio curved channel flow is by Kelleher et al. 

(Ref. 5]. These experiments show that, by traversing a hot- 
wire in a spanwise direction at several different radii, the 
secondary flow in the channel is periodic and undergoes a 
phase change as the probe is moved from the concave to the 
convex side of the channel. The range of Dean numbers for 
this experiment and accompanying flow visualization is 79 to 
PSK 

Using the same channel as in this study, Siedband 
(Ref. 6] and Niver [Ref. 7] used flow visualization to study 
Dean vortex structure and development. Siedband used a 
smoke wire at a Dean number range of 167-461 and at one 


streamwise location. Niver photographed Dean vortices 


using smoke at several different streamwise locations and a 
Dean number range of 40-220. The same 40 to 1 aspect ratio 
and radius ratio of 0.979 were employed in these studies as 
was used by Kelleher et al. (Ref. 5]. 

Finlay et al. [(Ref.1, 8] used the three-dimensional, 
incompressible Navier-Stokes equations to numerically 
Simulate curved channel flow with a geometry also similar to 
the Kelleher experiment. These simulations determined that 
the axisymmetric Dean vortices develop at least two 
different types of streamwise travelling waves. At a Dean 
number of 63.8 undulating Dean vortex flow consists of 
developing vortices and long streamwise wave lengths. At 
Dean numbers of 78.5 and 98.2, twisting Dean vortex flow is 
evident. This is characterized by a much shorter streamwise 
wave length and higher unsteadiness. Ligrani and Niver 
(Ref. 9] supported Finlay’s findings and concluded, through 
flow visualization, that Dean vortices develop such that 
they increase in radial height as the Dean number changes 
from 40 to 100. In this range, the vortices show some 
Spanwise and radial unsteadiness but are generally 
symmetric. Above a Dean number of about 100, vortices start 
to twist and are non-symmetric. Spanwise and radial 
oscillations become more pronounced and vortex flow becomes 
more unsteady. In all cases, the state of the vortex pairs 
1s dependent upon the streamwise location as well as the 


flow rate. 


B. OBJECTIVES 
The objectives of this thesis are; 


1. To measure the three mean velocity components of Dean 
vortices over a range of Dean numbers from 40 to 220. 


2. Using the results from objective 1, determine the 
three components of mean vorticity over the same 
Dean number range. 

3. To measure the frequencies of different unsteady 
events through the use of spectra of subminiature 
and conventional hot-wire voltage fluctuations. 


4. To investigate flow characteristics using flow 
visualization techniques. 


5. To compare qualitative and quantitative results 
with those given by Ligrani and Niver [Ref. 9], 
Kelleher et al. [Ref. 5] and Finlay et al. 

(Ref. 1, 8]. 

C. ORGANIZATION 

Subsequent to this introduction, Chapter 2 discusses the 
curved channel facility used in this study. Chapter 3 
explains, in detail, the five hole pressure probe used in 
the experiments to include its calibration. Chapter 4 
describes spatial resolution and downwash velocity 
corrections developed and adapted for five hole probe 
measurements. Chapter 5 shows the results of channel and 
probe qualification tests including qualification of 
pressure and velocity corrections. Chapter 6 describes 
subminiature and conventional hot-wire probe measurement 
equipment and procedures. Chapter 7 describes flow 


visualization equipment and techniques. Chapter 8 presents 


the results of this study and comparison with previous 


studies. Chapter 9 gives a summary and conclusions. 


II. CURVED CHANNEL 


A. CHANNEL 

Details on the design and construction of the Naval 
Postgraduate School Curved Channel Research Facility can be 
found in Siedband [Ref. 6], Niver [Ref. 7:pp. 10-13] and 
Ligrani and Niver [{Ref. 9]. A schematic [Ref. 7] and a 
photograph [Ref. 6] are shown in Figures 3 and 4. For 
completeness, a discussion on the channel and recent 
modifications are provided below. 

The channel is an open circuit suction facility designed 
for low speed transition studies. The inlet section 
consists of an aluminum honeycomb and three screens which 
reduce the spatial non-uniformities in the flow. The inlet 
is followed by a 20:1 contraction ratio nozzle. Flow is 
accelerated through the nozzle into a 2.44 m long straight 
duct. This entry length allows for fully developed channel 
flow entering the curved section for all Dean numbers 
studied. The flow then enters a curved test section with a 
convex wall radius of 59.69 cm and concave wall radius of 
60.96 cm. After the flow completes a 180 degree turn it 
enters another 2.44 m straight duct and goes through another 


honeycomb, screens and a diffuser into an outlet plenun. 


A 5.08 cm diameter pipe connects the outlet plenum to a 
second plenum which is kept at a low pressure by a 1/3 
horsepower blower. A globe valve for flow throttling and an 
orifice plate for mass flow measurement are located in the 
connecting pipe. The upper and lower walls of the channel 
are made of seamless sheets of polycarbonate (Lexan). The 
side walls, longitudinal supports and cross beam supports 
maintain a cross section of 1.27 cm by 50.8 cm within + 
0.015 cm for channel depth, and + 0.05 cm for the channel 
width. [Ref. 10] 

In order to isolate the channel from the blower and 
laboratory floor vibrations, rubber mounts were added 
between the inlet section and the plenum and their 
respective frames. Rubber vibration mounts were also used 
to mount the blower and blower frame. 

A slot was cut in the curved section of the channel to 
allow pressure measurements. The slot is 0.32 cm wide and 
7.62 cm long and is located in the convex wall of the 
channel. It is aligned in the spanwise direction 5.08 cm 
off of the centerline of the channel at the 120 degree 
streamwise location. Since cutting the slot at this 
location caused stress relief in the convex wall, care was 
taken to construct a support block which could eliminate any 
non-uniformities in the channel wall and provide an air 
tight seal through which a pressure probe could be inserted. 
A Lexan support block, (15.24 cm x 6.35 m X 2.22 cm) was 
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fabricated for this purpose. While the support block was 
being epoxied to the convex wall, the wall in the vicinity 
of the slot was pulled inward, toward the center of the 
curvature radius. This eliminated any bowing of the convex 
wall around the slot area. Foam strips were adhered to 
brass inserts which were placed just inside the slot. The 
brass inserts and compressed foam were positioned so that a 
probe could be inserted into the channel and move freely in 
the spanwise and radial directions without any air leakage. 
A schematic of the Lexan support block and the foam sealing 


device is shown in Figure 5. 


B. DEAN NUMBER DETERMINATION 

The mass flow rate and Dean number were calculated by 
measuring the pressure drop across the 1.5 inch ASME 
Orifice. The orifice pressure drop is measured using a 
Celesco model LCVR differential pressure transducer and a 
Celesco DC10D carrier demodulator. The signal is collected 
and processed using a Hewlett-Packard (HP) 3498A extender, 
HP3497A data acquisition system and a HP300 model 9836S 
computer. Additional details on the transducers, DAS and 
computer are found in the next chapter. 

Once the pressure drop is determined, a first guess at 
the Dean number is made by interpolating pressure drop vs. 


Dean number data from Niver [Ref. 7:p. 33]. Density is 


calculated using ambient temperature and pressure readings 
and the perfect gas relationship. The expansion 
coefficient, Y, is calculated from Holman and Gajda 

(Ref. 11}. Using the first guess Dean number to give an 
initial mass flow rate, the pipe Reynolds number, Re, is 


calculated as follows; 
Rep = ma,/pv A, 

The flow coefficient, K, is then interpolated from the ASME 
Tables [Ref. 12]. Mass flow rate is updated using the 
equation; 

m = KA,-Y¥ /2pAPoy 
This procedure is continued iteratively until the mass flow 
rate converges to within one percent of the previous value. 
Finally, the Dean number is calculated using; 


De = (m/p Ach) x (d/vemex 1 a7 


III. PRESSURE PROBE APPARATUS AND CALIBRATION 


A detailed investigation on the fabrication, use and 
calibration of five hole pressure probes is presented by 
Ligrani et al. [{Ref. 10]. te Seer ide a complete explanation 
of experimental apparatus and methods used for this thesis, 


portions of that article are presented in this chapter. 


A. BACKGROUND 

In the present investigation, a miniature five hole 
pressure probe having a probe tip diameter of 1.22 mm was 
developed to meet requirements of spatial resolution and 
minimal flow blockage. The probe was developed to measure 
three components of mean velocity at individual locations 
within the curved channel. In the working section of the 
channel, the cross-section is 12.7 mm x 508.0 mm. With this 
limited access for instrumentation, insertion and use of any 
existing type of triple-sensor hot-wire probe was not 
considered. The large meaSuring volumes and high cost of 
laser-doppler velocimetry systems made them an impractical 
alternative. 

Procedures for the calibration and application of 
ordinary-sized five hole pressure probes are described by 
Treaster and Yocum [Ref. 13]. Here, calibration 


coefficients for the yaw angle, pitch angle, and total 
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pressure are independent of Reynolds number based on probe 
tip diameters for Reynolds numbers ranging from 2000 to 
7000. The static pressure coefficient shows small Reynolds 
number dependence which, if unaccounted for, results in 
errors in magnitudes of measured velocity. Five hole probes 
similar in size to the ones used by Treaster and Yocum are 
used extensively throughout the fluid dynamics community and 
many results obtained from such probes are available in the 
literature. Treaster and Houtz [Ref. 14] describe the 
fabrication and calibration of miniature five hole pressure 
probes with tip diameters ranging from 1.07 mm to 1.68 mn. 
Calibration coefficients are defined in the same manner as 
given by Treaster and Yocum and also show similar dependence 
on Reynolds number and wall proximity. In a companion 
paper, Welz [Ref. 15] describes a double interpolation 
procedure to obtain velocity direction and magnitude from 


calibration data and measured pressures. 


B. PROBE DESCRIPTION 

Two schematics and a photograph of the miniature five 
hole pressure probe are shown in Figures 6, 7 and 8 
respectively. The diameter of the probe tip is 1.22 mn, 
where the central tube is surrounded by four other tubes, 
where the end of each is tapered at a 45 degree angle with 
respect to the mouth of the central tube. Each of the five 


hole tubes is made from stainless steel hypodermic tubing 


IBID 


with inner and outer diameters of .203 mm and .406 mn, 
respectively. The probe was constructed by first spot 
welding together these five tubes. The five were then bent 
in a 90 degree turn with a radius of 1.52 mm to avoid any 
tube crimping. Silver solder was then used to join the five 
tubes together continuously along their length. Finally, 
the portions of the probe tip were tapered to appropriate 
angles using a surface grinder. This approach is different 
from the one described by Treaster and Houtz wherein probe 
tips are ground and sanded into appropriate shapes after 
solder is used to completely encase all tubes. 

Figure 6 shows that the present probe tip extends 5.08 
mm (4 tip diameters) beyond the bend in the stem to minimize 
the effects of stem flow blockage at the measuring location. 
Similarly, the tip is 22.9 mm (19 tip diameters) from the 
9.53 mm diameter cylinder. Within this cylinder, each tube 
is connected to a 1.60 mm diameter tube using a small 
manifold for each tube connection. The mounting cylinder 
was sized for easy handling and mounting. On top of the 
cylinder is a second set of manifolds with five male 
connectors for attachment to 6.35 mm outer diameter/4.67 mm 
inner diameter polyflow tubing. 

The miniature five hole probe was designed to minimize 
the adverse effects of flow blockage and spatial resolution 


when used for measurements in the curved channel. When in 


12 


place for measurements at the centerline of the channel, the 
mounting cylinder is not exposed to the flow, and the probe 
produces 4.8 percent blockage over a 12.7 mm x 12.7 mm area. 
The tip is conical rather than a prism type to minimize 
near-flow blockage effects which are particularly severe for 


prism type probes. 


C. EXPERIMENTAL DETAILS 

1. Probe and Tubing 

The tubing used to connect the probe to transducers 

influences the pressure detected at the transducer. In 
particular, low amplitude pressure oscillations may result 
when there is a large volume of air in the tubes relative to 
that in the probe. With zero differential pressure applied 
across the transducer, these oscillations have peak to peak 
pressure amplitudes about .0003 cm of water. The period of 
oscillation is approximately one minute. These variations 
are believed to result due to compressibility of the gas and 
the constriction provided to the air as it passes through 
the .20 mm inner diameter holes in the probe tip from the 
larger diameter tubing. The motion of the air and the 
constriction result in pressure waves which seem to initiate 
from small disturbances to the tubing. As the volume of air 
in the tubes increases, the pressure detected at the 
transducer becomes more unsteady and longer decay times are 


required before the measurement reaches a steady state. 


1b) 


Ordinarily, to overcome such effects, tubing is fixed in 
place so that it may experience minimal vibration. Another 
alternative, used in the present study, is to minimize the 
volume of air in tubes. Flexible plastic tubing with 

3.16 mm outer diameter and 2.66 mm inner diameter 
(commercially known as 1/8 inch Polyflow), less than 20 cm 
long, is employed for this purpose along with plastic 
adapters to connect the tubes to the male connections of the 
probe manifold. For differential pressures less than 0.2 cm 
of water, the pressure unsteadiness is then reduced 
considerably. 


2. Transducers, Carrier Demodulators and 
Signal Processing 


Two different sets of transducers are used to 
convert pressure to an electrical signal. For differential 
pressures less than .25 cm of water, Validyne model DP103-06 
variable reluctance differential transducers are employed. 
For higher pressures, Celesco model LCVR variable reluctance 
transducers are used. A separate transducer is used to 
measure pressure relative to the atmosphere for each port of 
the five hole probe. 

The five Validyne transducers each have a full scale 
range of .25 cm of water differential pressure, and are used 
to measure pressures as low as .013 cm of water or five 
percent of the full scale range. These devices are quite 


sensitive, requiring extra care as they are used, 
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particularly as tubing from the probe is connected or 
disconnected. This is done very slowly while monitoring 
transducer output to avoid pressure overload. Each 
transducer is supported in a vertical position by a 
plexiglass mount slightly above the probe. This mount 
minimizes gravitational effects, vibration and orients 
transducers to minimize the influence of magnetic fields 
from two nearby stepping motors. Shielded cables connect 
the transducers to Celesco CD10D carrier demodulator units. 
The shields extend to the transducer cases to provide a 
ground connection to the demodulators and significantly 
reduce electronic noise. Further reductions are achieved by 
mounting the demodulators on a thick metal plate to further 
isolate them from the stepping motors. The carrier 
demodulator units are essentially AC bridge circuits which 
produce a DC voltage output proportional to the pressure 
Signal from the transducers. 

The signal to noise ratio of the Validyne 
transducer/carrier demodulator system ranges from 100 to 
5000. The amount of temperature drift for zero differential 
pressure is typically less than .01 volts over an eight Hone 
period. This drift was recorded to be as high as about .1 
volt, but only for one transducer. This is still quite low 
considering that the maximum voltage produced from a 


demodulator is 10 volts. Demodulators are adjusted to 
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maximize output voltage in order to minimize the adverse 
effects of zero drift and to maximize the signal to noise 
evel ALS) « 

The Celesco model LCVR variable reluctance 
differential pressure transducers are supported by the short 
tubes extending from the neuen of the probe. The full scale 
range of each transducer is 2 cm of water. Celesco carrier 
demodulators are again employed to convert electrical output 
to DC voltage. When compared to the Validyne transducers, 
Celesco devices have higher noise levels and for 
differential pressures less than .2 cm of water, more 
Significant calibration non-linearities. 

Each transducer/demodulator system is calibrated by 
first adjusting the zero so that the output with no applied 
pressure is less than .01 volts. This small offset voltage 
is determined from the average of 50 samples taken by the 
data acquisition system, recorded and then accounted for in 
subsequent Signal processing. Next, a constant pressure, 
close to the value expected during measurement, is applied 
Simultaneously to the unused orifice of the transducer and a 
Meriam Model 40GD10 WM 1.27 cm horizontal manometer, which 
can be read to an accuracy of 0:°5005 cm enhevacen ee. 
calibration coefficient is then determined to be the ratio 
of the manometer pressure and the difference between the 
average voltage and offset voltage. This procedure is 
repeated three times to give three calibration coefficients 
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which are then averaged to reduce experimental uncertainty. 
A typical value for Celesco transducers is 0.15 cm of water 
differential pressure per volt. Calibration coefficients 
for Validyne transducers are typically about .015 cm of 
water per volt. 

3. Computer and Software 

Each transducer/demodulator system is connected to a 
Hewlett-Packard (HP) 3498A extender which is controlled by a 
HP3497A data acquisition/control unit. The HP3497A is 
equipped with analog multiplexing and a digital voltmeter 
with l»V sensitivity. The 3497A is controlled by a HP 
Series 300, Model 9836S computer, which also processes the 
data. The HP9836S is equipped with a MC68000, 8MHz 16/32- 
bit processor, dual 5% inch floppy disk drives, and a memory 
of 1M bytes. 

When data is acquisitioned from the probe, voltage 
samples are recorded sequentially from the five channels 
until 50 samples have accumulated for each channel. The 
mean and standard deviation are subsequently calculated for 
each channel. 100 additional samples are then recorded for 
each channel provided each is within three standard 
deviations of the current mean value. The overall sampling 


rate is about 5Hz. 
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D. PROBE CALIBRATION 

1. Wind Tunnel 

The wind tunnel used for probe calibration is 

located in the laboratories of the Department of Mechanical 
Engineering of the Naval Postgraduate School. It is an open 
circuit blower tunnel which provides uniform flow at the 
nozzle exit with a free stream turbulence intensity of 
approximately 0.1 percent at 30 m/s. Air enters the fan 
inlet from the surrounding air through a coarse filter. A 
1.6 mm clearance between the fan and diffuser isolates 
vibrations from the fan to the wind tunnel body. The 
diffuser contains a second fine filter to remove small 
particles from the air as well as four baffle vanes to 
reduce noise and minimize the likelihood of flow 
separations. This diffuser is then followed by a header 
containing a honeycomb and three screens to reduce spatial 
nonuniformities. A 16:1 contraction ratio nozzle leads to 
the 3.05 m long, 0.610 m wide test section. The height of 
the top wall is adjustable. For this work, a zero pressure 
gradient is maintained to within 0.013 cm of water 
differential pressure along the length of the test section. 


The exit of the tunnel vents to the laboratory. 
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2. Calibration sled 

To set probe angles during calibration, the probe is 
positioned in a manually operated traversing unit 
manufactured by United Sensor and Control Corp. The unit 
has been modified to include a protractor to facilitate 
positioning the probe to an accuracy of 0.5 degrees over a 
range of yaw angles from -40 degrees to +40 degrees. The 
traversing device is mounted on a horizontal sled which fits 
on top of the wind tunnel side walls. The mounting sled can 
be rotated through a range of pitch angles from -45 to +45 
degrees to an accuracy of 0.5 degrees. The pressure 
transducers rotate with the probe and are sensitive to the 
direction of gravity and thus the pitch angle. To account 
for diaphragm deflection from gravity as the transducers are 
rotated, a constant multiplied by the sine of the pitch 
angle is added to the pressure signal. 

Figure 9 shows a photograph of the calibration sled 
mounted on the top of the wind tunnel with probe, traversing 
device, protractor and transducers. 

3. Calibration Procedure 

The probe is calibrated in the uniform free stream 
of the wind tunnel. The right handed coordinate system 
shown in Figure 6 is used. Positive yaw angles exist when 
the pressure at port two is greater than that at port three; 
pitch is considered positive when the pressure at port four 


1s greater than that at port five. 
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The probe is first aligned in the wind tunnel. Zero 
yaw is defined to be that angle where the pressures at ports 
two and three are equal. The calibration sled is 
constructed so that at zero pitch angle, the probe mounting 
cylinder is vertical within 0.5 degrees. Here, the 
difference in the pressures at ports four and five produces 
a nonzero pitch coefficient, which 1s also used for zero 
pitch when measurements are made. 

To obtain calibration data, the probe is manually 
fixed at a known yaw angle and then rotated through the 
pitch plane. Similar procedures are employed by Treaster 
and Yocum [{Ref. 13] and Treaster and Houtz [Ref. 14}. A 
total of 25 angle combinations are used, five pitch angles 
at five different yaw angles. In addition to the pressures 
at the probe’s five holes, the total and static pressures 
are required. The total pressure is measured with a Kiel- 
type total pressure probe; the static pressure is measured 
at the same streamwise position with a pressure tap in the 
Side wall of the tunnel. To account for small changes in 
the free stream velocity, the static and total pressures are 
updated periodically. Typically, they vary less than 


0.0005 cm of water during a three hour calibration session. 
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For each probe position during calibration, yaw, 
pitch, total, and total minus static pressure coefficients 
are calculated using; 

2 195) 02 P) 


Soe = Wy 2 PSO = P) 


eee re - AP) = P) 
and 

Cots = (Pe - Pg)/(Py ~ P) 
respectively. Cop and Coy have the same form as the 
coefficients employed by Treaster and Yocum [{Ref. 13]. New 
expressions for Cot and Cots are used to produce 


coefficients which do not approach zero as the yaw angle and 


pitch angle approach zero. 


E. CALIBRATION RESULTS 
Calibration coefficients are presented in Figures 10-14 
for a mean speed of 1.9 m/s, which corresponds to a Reynolds 


eng oA yal 


number based on probe tip diameter Reprobe py 


Figure 10 increases with yaw angle and is close to zero at 
zero yaw angle. Some slight dependence on pitch angle is 
evident. Cop Variation with yaw and pitch in Figure 11 
shows similar behavior, except it is about equal to -.15 at 
zero pitch angle for all values of the yaw angle tested. In 


Figures 12-14 Cote Cots and Cos show significant dependence 


on yaw angle at a given pitch angle. Here, 


Pad 


C = 2 eee) Eee P. The Cot and Cots curves are somewhat 
non-symmetric about zero pitch angle indicating a slightly 
non-symmetric probe tip. 

The dependence of the yaw coefficient on probe Reynolds 
number is illustrated by results presented in Figures 15-18. 
These results correspond to mean speeds ranging from 1.1 to 
5.0 m/s and Reynolds numbers ranging from 86 to 391. Except 
for some differences for Re~86 due to scatter resulting from 
measurements of extremely low pressure levels, the yaw 
coefficient calibration curves are repeatable with minimal 
scatter and are independent of Reynolds number. This 
includes the value of Coy at zero yaw angle and zero pitch 
angle. Other coefficients show similar behavior, in 
contrast to the results of Treaster and Yocum wherein the 
static pressure coefficient shows significant Re€probe 
dependence. The forms of the presently used coefficients 
were chosen to avoid this problem. 

Calibration results and data reduction procedures were 
checked by placing the probe in the wind tunnel at known yaw 
and pitch angles and then determining these angles using the 
measurement procedures. At a mean speed of 2.7 m/s, both 
measured angles were generally within .3 degrees of preset 


angles. Differences of about .5 degrees occurred only when 


the yaw angle was greater than 15 degrees. 
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Based on a 95 percent confidence level, the 


uncertainties of Coy and Cop are both estimated to be about 


+ .02. The uncertainty of C iSs-typlcally + .042) The 


pts 
uncertainties of the slopes of the yaw and pitch calibration 
curves are approximately + 4.4 percent. A sample 


uncertainty analysis is given in Appendix B. 
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IV. MEASUREMENT CORRECTION PROCEDURES 


Adjustments must be made to account for pressure and 
velocity gradients measured in a three-dimensional flow 
field. Two of the most important adjustments result from 
Spatial resolution and induced downwash velocity 
corrections. Ligrani et al. [Ref. 16] present a detailed 
account of the literature reviewed and methods used to 
employ the corrections used in this study. The results of 
this study are presented in this chapter as they apply to 


curved channel measurements. 


A. SPATIAL RESOLUTION CORRECTION 

Accounting for spatial resolution of multi-hole pressure 
probes is necessary since different pressure ports are at 
different locations. The simplest accounting corrects 
pressures so that all appear to be measured at the same 
physical location. For five hole probes, if this common 
point is chosen to be the location of the central hole, then 
four pressures require adjustment. 

Figure 19 illustrates how the spatial resolution scheme 
is applied. Here, locations of the five pressure ports are 
shown as the probe is moved to different z locations at a 
constant value of y. For this case, to simplify discussion, 


corrections are applied only to ports two and three. In 


24 


most real applications, the scheme is also used 
Simultaneously for ports four and five, provided probe data 
are available for at least three locations in different 
planes normal to the wall. Spanwise gradients in total 
pressure appear as differences in pressures measured at 
ports two and three. Thus, for the situation of Figure 19, 
pressures from ports one, four and five are not affected by 
spanwise gradients. Only P5 and P3 require adjustment. 

The correction is implemented after pressures from each 
of the five holes are recorded at measurement locations 
k=1, 2, 3, ... n. As a first step, cubic spline curves are 
fit to pressures measured at port two for probe positions 
k=1, 2, 3... n. The exact spanwise position of this port 
is used and not the location of the center hole. In Figure 
19, these locations of P5 measurement are z, + d, where d is 
the distance between the center hole (port one) and the 
center position of each adjacent hole. A similar procedure 
is used for P, measured from port three, except that its 
locations are z} - dad. One smooth curve is constructed for 
each port, each of which is then used to calculate the value 
Same, and P, at the location of port one. In Figure 19, 
these pressures are designated P.,, and P3 ok! and evaluated 
at locations 2; and fully corrected to account for finite 
spatial resolution. They may then be used to calculate the 
yaw coefficient for determination of the yaw angle at 
location 2}. If pressures from ports 4 and 5 had been 
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corrected to give P4k and Ps x! then these pressures could 
then be used to calculate pitch coefficients for pitch angle 
determination at 2,. 

With this approach there are no adjustable coefficients. 
Required are y and z locations of port one and d, the 
distance between port one and the center of the other probe 
pressure ports. For the custom-made miniature probe, 
d= .041 cm. The scheme assumes that pressures vary 
smoothly without significant discontinuities between 
pressure port locations and measurement locations. At the 
end points of any traverse, satisfactory results may be 
obtained by extrapolation. For optimal accuracy, in the 
present study, endpoints are discarded after interpolation 
for their neighbors is completed. 

This scheme is physically more realistic than others 
which employ streamwise velocity gradients as a basis to 
account for spatial resolution. Finite spatial resolution 
is caused by variations of pressure, and thus, pressures 


must be corrected and not velocities or velocity gradients. 


B. INDUCED DOWNWASH VELOCITY CORRECTION 

The downwash velocity is part of induced secondary flow 
which occurs as a result of the blockage of blunt bodies 
within transverse gradients of mean streamwise velocity. 


Young and Maas [Ref. 17] show that velocities measured in 
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boundary layers using flat-faced pitot probes correspond to 
positions farther from the wall than where the measurements 
are actually made. Livesey [Ref. 18] obtains similar 
results for flat-faced probes, but finds that apparent 
streamline displacement depends strongly on probe shape. 
Lighthill [Ref. 19] describes a theoretical model for flow 
around a sphere in a uniform velocity gradient. He finds 
that the displacement of the stagnation streamline is caused 
by downwash velocity near the sphere with flow from regions 
of high streamwise velocity to regions of lower streamwise 
velocity. He also confirms earlier work showing that 
downwash velocity can be quite significant and depends 
strongly on the shape of the probe. For multiple hole 
probes that are used to determine flow direction as well as 
speed, the induced downwash velocity, if uncorrected, causes 
larger errors than displacement of the stagnation 
streamline. 

After pressures and pressure coefficients are corrected 
for effects of finite spatial resolution, the velocity 
distribution is determined as described in Chapter 5. With 
this velocity field in the computer memory, adjustment for 
downwash velocity effects can proceed. 

The present correction for induced downwash velocity is 
based on the work of Lighthill (Ref. 19]. He shows that 
D/(a’aU/ds) is of the order of one near a sphere. Here, a’ 


is the sphere radius, and D is the induced downwash velocity 


27 


in the s direction, which is normal to the direction ome 
adU/ds is the gradient of the U component of velocity with 
respect to the s-direction. For probe tips, downwash 
velocities in the normal and spanwise directions may be 
estimated using; 


du,, du 


D., = --- a’E and D, = --- a’/E 
z, 
Yay Y dz : 


where a’ is the probe tip radius and dimensionless quantities 
Ey and E, are coefficients which depend on probe geometry. 
These may be different from each other, especially for non- 
axisymmetric probe tips. In the present study, however, 
these two coefficients for the normal and spanwise 
directions are assumed to be the same and equal to E. 


Normal and spanwise velocity components U,, and U, corrected 


x 
for the effects of downwash velocity are then given by; 


au au 
U,, = U ti > 4 and U, = U 5 — oa 
u Zz Zi 
yy ¥8 ay dz 
where A=a’E. In the above two equations, Uvu and Uri 


represent velocities uncorrected for the downwash effect. 
qU,/dy and dU,/dz are determined by fitting U, with cubic 
spline curves in the y and z directions. The cubic splines 
are differentiated analytically and the values of the 


derivatives at every measurement point are then calculated. 
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In the present application, the final velocity 
components are corrected by numerically fitting the cubic 
splines and differentiating them. The corrected velocities 
are stored to be used for later plotting and vorticity 


computations. 
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V. CURVED CHANNEL MEASUREMENT AND QUALIFICATION 


This chapter describes measurement procedures and 
channel qualification tests which are presented by Ligrani 
et al. (Ref. 10, 17]. Portions of these articles are 


presented here as they apply to this study. 


A. MEASUREMENT PROCEDURES AND APPARATUS 
1. Traversing Mechanism 

When meaSurements are made in the channel, the probe 
is oriented using a mount attached to a travelling block 
which is part of an automated two-dimensional traversing 
mechanism. The set-up is shown in Figure 20. The 
travelling block moves in the spanwise and radial 
directions. These parts are driven by 20-thread per inch 
drive screws and guided by ground steel, case-hardened 
Support shafts. The drive screws are rotated by Superior 
Electric type M092-FD310 stepping motors. These motors are 
controlled by a MITAS Motion Controller, which directs the 
movement of the probe according to instructions from the 


HP9836S computer by means of an RS232C data interface. 
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The mount used to orient the probe is also evident 
in Figure 20. This mount consists of two collars. The 
first may be rotated relative to the radially mounted 
travelling block to allow the pitch angle of the probe to be 
set relative to the channel flow. A second collar is 
attached to the probe by means of a set screw, which allows 
the yaw angle of the probe to be adjusted. Changes to the 
yaw angle as small as +.1 degrees and to the pitch angle as 
small as +.5 degrees are possible with this mechanism. 

2. Measurement Procedure 

The probe is positioned and aligned before 
measurements are made. The probe is first inserted into the 
channel at a spanwise location 5.08 cm from the centerline 
so that part just touches the concave wall. Then it is 
visually aligned with the flow and subsequently moved to the 
channel centerline using the MITAS controller. Small 
adjustments to the probe yaw angle and pitch angle are now 
made using the probe mount to obtain a zero yaw coefficient 
and a small pitch coefficient. A bulk flow velocity of 
0.3 m/s is used for alignment since flow here is in the 
streamwise direction and minimal secondary flows are 
present. The low speed also means that the pressure 
differences between the ports are small, hence there is some 
scatter in the pitch and yaw coefficients. The final 
orientation of the probe gives a yaw coefficient of 


0.0 +0.01 and a pitch coefficient of -.60 + .01. This is 


om 


equivalent to +4.5 degrees offset accounting for the pitch 
coefficient of -.15 obtained during calibration 

alignment. More accurate pitch alignment forces the probe 
against the edge of the slot, hence the pitch offset is 
corrected during data reduction. No offset is applied or 
needed for the yaw angle. 

In a typical spanwise/radial plane traverse, the 
probe moves in the radial direction from the concave wall to 
the convex wall as measurements are made. To begin 
measurements at a new spanwise location, the probe first 
moves 0.254 cm in the spanwise direction, then in the 
reverse span direction 0.127 cm, and finally toward the 
convex wall until it is .254 cm away. A new radial 
measuring traverse then begins. This procedure is used for 
positioning to reduce the adverse effects of backlash and 
give repeatable data regardless of the direction of spanwise 
probe motion. 

When measurements are made at a particular location, 
150 samples of each pressure are averaged. This is done by 
Sampling pressure signals from probe ports in sequence so 
that a longer averaging time is provided for each channel, 
thereby minimizing the influence of short duration 
unsteadiness in the channel. -For each probe location, about 
two minutes of acquisition and processing time is required. 


A traverse over an entire plane requires about 
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eight hours. The entire process is computer controlled and 
generally performed at night when there is minimal 
unsteadiness at the channel inlet resulting from laboratory 
disturbances. 
3. Data Reduction 

Measure pressures from each probe port are stored by 
the computer until an experiment 1s completed. These data 
are then corrected to account for the effect of finite 
spatial resolution using procedures described in Chapter 4. 

The next step is to determine the values of C and 


Py 


Cop for each location using measured pressures. A four step 
procedure is used to determine the yaw angle and the pitch 
angle. First, a preliminary pitch angle p’ is calculated 
using; 

0 Oo (SG aeach Mey = jee: 
This equation provides a reasonable fit to pitch angle 
calibration data such as that shown in Figure 11 for a mean 
velocity of 1.9 m/sec. Here, Por¢e=4-5 degrees, the pitch 
angle offset resulting from the probe mount. Referring to 
yaw angle calibration data in Figure 10, the second step 
requires the determiniation of a preliminary yaw angle y’ by 
interpolation using p’ and the measured value of Choy: 19! 
the third step, pitch angle is determined using the measured 
value of Cop and y’, again interpolating from results such 
as the ones presented in Figure 11. The final pitch angle p 


is then determined by subtracting Poff: In the 
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fourth step, the final yaw angle y is determined from Coy 
and p by interpolating results such as those in Figure 10. 
Now that the yaw angle and pitch angle are known, 


E is computed using a double linear interpolation 


pts 


procedure. A plot of C 


pts versus yaw angle and pitch angle 


at a mean speed of 1.9 m/s is shown in Figure 13. Knowing 
Cots: the total velocity magnitude U is determined using; 
U = (2 Copy (Py - P)/p)? 

where p is the fluid density. This approach gives less 
uncertainty for U than one which employs two separate 
coefficients for the total pressure and static pressure 
(l.e. Treaster and Yocum, [Ref. 14]). The total pressure 
relative to atmospheric pressure P,,, at each measurement 
location then follows using; 

Pamb ~ Pt = (Pamb ~ P) ~ Cpt (Pi ~ P) 
Cot 1s obtained using a double linear interpolation scheme 
from results such as the ones given in Figure 12. A 
coefficient for static pressure C is not used for data 


ps 


reduction. The three velocity components are subsequently 


given by; 
Ug = U cos(yaw) cos(pitch) 
U, = U sin(yaw) 
U, = U cos(yaw) sin(pitch) 


In the final phase of data reduction, corrections 
are made to account for the downwash velocity correction as 


described in Chapter 4. 
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B. CHANNEL QUALIFICATION TESTS 

In order to qualify the performance of the curved 
channel, a number of tests were conducted. Many of these 
tests are explained in detail by Ligrani and Niver [Ref. 9] 
and are summarized here. 

Turbulence intensity levels were measured at the nozzle 
exit using a horizontal hot-wire probe. The probe was 
mounted normal to the flow direction and operated in the 
constant temperature mode. At a Dean number of 290 (local 
mean velocity of 3.4 m/sec) a turbulence intensity of 0.0006 
was measured. At a Dean number of 450 (local mean velocity 
of 5.2 m/sec) a turbulence intensity of 0.001 was measured. 

A miniature Kiel probe was used to ensure that the 
geometry of the channel at the measurement slot was uniform. 
The Kiel probe provides an accurate measure of total 
pressure when flow vectors are misaligned as much as 40 
degrees from the probe tip. The diameters of the shield and 
pressure port of this probe are 1.59 mm and 0.4 mm. The 
same traversing mechanism and procedure is used for Kiel 
probe measurements as was described for the five hole probe. 
The Celesco 2 cm transducer was used to obtain measurements 
with the Kiel probe. Results of total pressure measurements 
relative to ambient pressure for a Dean number of 54 are 


shown in Figure 2la. 
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Several different tests were made to determine the 
effects of perturbations on the flow field. BlockKages of as 
high as 50 percent of the inlet nozzle did not appear to 
disturb the measurements taken at the 120 degree location or 
cause flow visualization results at any point in the curved 
section to be different. The effect of slight non- 
symmetries or discontinuities were checked by temporarily 
deforming the convex and concave walls of the channel or 
placing a strip of tape near the nozzle exit. No changes to 
the flow were noted for Dean numbers less than 200. It was 
found that leaks of any size, particularly upstream of, and 
near the location of measurement points would have a 
Significant effect on measured results. Unsteadiness forced 
on the flow by moving objects or laboratory personnel around 
the nozzle entrance also had an effect on measured results. 
To minimize the effect of unsteadiness on experiments, they 
were run at night when there was little activity in the 


laboratory spaces. 


C. FIVE-HOLE PROBE QUALIFICATION 

In order to qualify the behavior of the miniature five 
hole pressure probe, the total pressure distribution was 
measured and compared with measurements from the miniature 
Kiel probe. Total pressure measurements for a Dean number 
of 49 are shown in Figure 21b. Comparison with the Kiel 


probe measurement (Figure 21a) shows excellent agreement 
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providing evidence of accurate five hole probe alignment and 
measurement procedures. Comparison is made at this Dean 
number since minimal secondary flows are present and because 
gradients of mean streamwise velocity are very small in the 
spanwise direction. For the data presented in Figure 21a, 
40 spanwise locations and 8 vertical ones are used with a 
resolution of 0.127 cm. The measurement area begins 
0.127 cm from the concave wall, 5.08 cm from the channel 
centerline and extends to 0.254 cm from the convex wall, 
10.16 cm from the channel centerline. 
D. SPATIAL RESOLUTION AND DOWNWASH VELOCITY 

CORRECTION APPLICATIONS 

Figure 22 shows measurements of streamwise velocity 
after corrections are applied, plotted at three radial 
traverse locations and compared to Curved Channel Poiseuille 
Flow (CCPF) profiles. Equations for CCPF are given by 
Finlay et al., [{Ref. 1]. The measurements were taken at a 
Dean number of 36 and are generally in excellent agreement 
with the CCPF profile for this Dean number. Small 
differences may be attributed to scatter caused by the 
extreme difficulty in measuring the low pressure levels 
associated with these velocities. 

Another example of the application of pressure and 
velocity corrections is now presented. At a Dean number of 


100, secondary flow velocities are reasonably steady and 
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Spanwise periodic. The velocity field thus has pressure 
gradient in both spanwise and radial directions. 

Figure 23a shows contours of streamwise velocity. This 
distribution depends only weakly on the corrections. Figure 
23b shows secondary velocity vectors when no corrections are 
made. Velocity vectors point away from regions of high 
streamwise velocity to lower speed regions. Such behavior 
is present for radial and spanwise components and is 
physically unrealistic, especially near surfaces. Figure 
24a shows secondary flow vectors from the same flow after 
the pressures are corrected only for spatial resolution. 
These vectors are substantially different from those in 
Figure 23b, even though vectors near walls still point away 
from the regions of high streamwise velocity. The largest 
vector in Figure 23b is .14 m/s, whereas the largest vector 
in Figure 24a is .10 m/s. In Figure 24b, both spatial 
resolution and downwash corrections are applied. Here, the 
value of A employed in the induced downwash velocity 
corrections is .025 cm or .4a’. Velocity vectors follow 
expected behavior. Well defined vortices are evident and 
vectors show no consistent tendency to point into walls or 
away from regions of high streamwise velocity. Secondary 
flow vectors indicate upwash from the concave to the convex 
surfaces over the same areas where streamwise velocity 


deficits are present in Figure 23a. 
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VI. POWER SPECTRA APPARATUS AND PROCEDURES 


Spectra of voltage fluctuations from hot-wire 
anemometers which resulted from velocity fluctuations at the 
hot-wire probe tip are now presented. Measurements were 
made by inserting a subminiature or conventional horizontal 
hot-wire into the curved channel at the 110 degree 
streamwise location. These hot-wire probes were mounted 
through a plug in the channel side wall and aligned so that 
wires were sensitive to streamwise and spanwise flow 
components in the case of the horizontal wire, or to all 
three flow components in the case of the subminiature 
crossed-wire probe. The hot-wires were mounted on the 
centerline between the concave and convex walls so that the 
wires were 9.60 or 9.92 cm from the channel side wall. Flow 
visualization has shown that this is an adequate distance to 
avoid side wall effects. A photograph of the horizontal 


wire probe mounting is shown in Figure 25. 


A. PROBES AND MEASUREMENT CHAIN 
1. Hot-Wire Probes 
The primary device used for spectral measurements 
was the subminiature crossed-wire probe. This was chosen to 
take advantage of the excellent spatial resolution 


associated with this type probe. This hot-wire was 
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manufactured at NASA Ames Research Center using fabrication 
techniques discussed in Westphal et al. [Ref. 20] and 
Ligrani and Bradshaw [Ref. 21]. The crossed-wire probe is 
manufactured of 90% platinum / 10% rhodium wire with a 
nominal diameter of 0.625 microns and a length of about 200 
microns. | 

The horizontal hot-wire probe used is a DANTEC type 
55P05. This probe is fabricated of platinum plated tungsten 
and has a diameter and wire length of 5 microns and 1.25 mm 
respectively. The horizontal wire was used to determine 
velocity fluctuations averaged over a much wider portion of 
the channel depth than the subminiature wire and, in some 
cases, resulted in slightly different spectra than from the 
subminiature wire. 

Checks were made to assure that no vortex shedding 
from hot-wires existed to affect power spectra. 
H. Schlichting (Ref. 22, pp. 28-32] presents a relationship 
between the Strouhal number and Reynolds number. Knowing 
the Reynolds number, one can determine the frequency of 
vortex shedding off of the wire, if it exists. Even when 
the diameter of the large wire and the maximum velocity 
measured are considered, the wire diameter Reynolds number 
for the present study is less than one, which indicates that 
there is no vortex shedding off of the hot-wire filaments. 
AS an additional check, vortex shedding off of the prongs 
supporting the wire filament was investigated. For the 
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large hot-wire the maximum Reynolds number attained is 4.9. 
Based on these calculations, it was determined that there is 
no vortex shedding off of the prongs. 

2. Signal Measurement 

The continuous signal from either hot-wire was sent 
to a DANTEC 55M10 CTA Standard Bridge via a five meter 
coaxial cable. The bridge was operated in the constant 
temperature mode. The subminiature hot-wire sensor has much 
higher resistance than a standard hot-wire. Due to this, 
the DANTEC bridge was modified by adding a 2000 ohm metal- 
foil resistor at the location of the .09 ohm setting 
resistor. Recalibration of the bridge against an ohmmeter 
indicated that 98.90 ohms should be added to the bridge 
reading when the decade controller was switched to the .09 
ohm position. Adjustment of cable impedence was 
accomplished by setting the cable compensator to the 5 m 
position. For subminiature hot-wire measurements, the 
bridge gain was set at "1". For the conventional horizontal 
wire the bridge gain was set at "2". Subminiature wire and 
conventional wire sensors were used at an overheat ratio of 
1.2 and 1.5 respectively. 

A DANTEC 56N20 Signal Conditioner was employed to 
filter and amplify the hot-wire signal. The HF filter was 
set to 0.1 Hz in all cases to remove the DC signal. To 
prevent aliasing, the LF filter was set at 100 or 300 Hz, 


depending on the frequency range being investigated. Gain 
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settings were maximized depending on the type of wire and 
flow rate. Typically, a gain of 100 was used with the 
subminiature wire at all flow conditions. A gain of 500 was 
used with the conventional wire at Dean numbers less than 
120. Above this flow rate a gain of 100 was used. 

Output from the DANTEC Signal Conditioner was sent 
simultaneously to signal processor equipment and a PRECISION 
model 1476A 10KHz oscilliscope. The oscilliscope was used 
to help determine the frequency ranges that should be 
investigated. 

For the 1.25 mm wire, the time response of the 
sensor-bridge combination was adjusted by using the square 
wave test function of the bridge. Using an overheat ratio 
of 1.5, a response time of 12.5 y»wsec was achieved. 
Photographs of the probe measurement chain equipment are 


shown in Figures 26a and 26b. 


B. SIGNAL PROCESSING AND ANALYSIS 

Two different methods were used to obtain spectra of 
Signals from the hot wire sensor. The first of these was a 
Hewlett-Packard high speed data acquisition system with 
accompanying Fast Fourier Transform software. The second 
method employed was a Scientific Atlanta real time signal 
analyzer. The two methods wane used simultaneously to check 
software and to identify spectral variations from digital 


and electronic noise. 
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1. High Speed Digital Acquisition and Analysis 

A HP6944A Multiprogrammer mainframe unit was used 
for digital processing of the hot-wire time domain signal. 
This utilizes a four channel HP69755A FET high speed data 
acquisition card in conjunction with a HP69750A scanner 
control card. Each channel has a maximum sampling rate of 
25 KHz. One channel operation is used for our purposes. A 
HP69759A 500 KHz, 12 bit analog to digital converter 
digitizes and passes data via HP14717A buffered ADC cable 
and HP69791A 64K RAM to the HP9836S computer. Data is 
gathered and sent to the computer as blocks of 2048 integer 
data points. A complete sampling of data consists of 
Sixteen blocks of data. The acquisition of data, including 
acquisition rate, is controlled by the computer using 
parameters input by the user. The digitized time signal is 
stored on floppy disk for analysis at a later time. 

The time domain signal is transformed to frequency 
domain using a Fast Fourier Transform routine adapted from 
Brigham [Ref. 23]. The transform and acquisition software 
were written in BASIC 4.0 by Dr. B.A. Singer. The FFT reads 
and processes each of the sixteen blocks data separately and 
ensemble averages the sixteen blocks. The Hanning Window 
weighting parameter is used to treat the side lobes for each 
block of data. Before plotting, the amplitude of the signal 


is converted to normalized power as (Volts?/Hz) so that the 
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integral of the power over the frequency range is one. The 
digital processing equipment and transform were tested using 
a DYNASCAN model 3010 Function Generator. 
2. Real Time Signal Analysis 

For real time analysis of the hot-wire signal, a 
Scientific Atlanta (Spectral Dynamics Division) SD380Z 
Signal Analyzer was used. The SD380Z is a four channel, 
microprocessor based Signal analyzer. The analyzer was used 
in the spectrum mode of display with the trigger set for 
"free run". This allows data acquisition at the maximum 
rate depending on the frequency resolution used. For the 
800 line frequency resolution used, this equates to a sample 
size of 2048 bits of data sampled at 10 KHz. The Hanning 
Window FFT weighting parameter was chosen so that both 
methods of signal analysis could be compared. Similarly, 
data was plotted as Volts“/Hz vs. Hz. Blocks of transformed 
data were ensemble averaged over the maximum of 999 
ensembles. A HP7470 plotter was used with an IEEE 488 
interface to plot the spectra display. The advantage to 
using this signal analyzer is that it allows for much faster 
processing and display of spectral plots which allows a 
large number of frequency ranges to be investigated ina 
short time. Data is presented in Chapter 8 using both real 
time signal analysis and the HP/DAS-FFT system. It is 
evident from the plots (Figures 75-120), that the systems 


produce nearly identical results. 
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VII. FLOW VISUALIZATION APPARATUS AND PROCEDURES 


Flow visualization photographs were obtained to 
illustrate the qualitative behavior of Dean vortices. Both 
the radial-spanwise plane and the exterior of the channel 
(viewed from the outer, or concave side) were observed and 
photographed. Equipment and methods used are presented 
below. 

A. SMOKE GENERATOR 

Rather than using the wood-burning generator used by 
Niver [Ref. 7], a commercial fog/smoke generator was 
employed. This device was selected because it is safer, 
more compact and more convenient than the wood burning smoke 
generator. The commercial generator is a ROSCO model 1500 
Fog and Smoke Machine, which uses ROSCO Fog/Smoke Fluid, a 
mix of glycols of low molecular weight. The fluid is drawn 
into the smoke generator, heated to near its vaporization 
point and atomized by forcing the fluid through an orifice. 
The smoke is then pumped into an inlet plenum and injected 
into the entrance of the channel through 1.27 cm inside 
diameter tubes. The smoke plenum was designed to provide 


laminar jets of smoke at the entrance of the curved channel. 
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The best method for using the smoke generator for flow 
visualization was to inject a puff of smoke into the channel 
entrance and then turn the smoke generator off. After a 
period of time a steady stream of smoke could be viewed at 
the location to be photographed. Generally, pictures could 
be taken over about a ten second period. The amount of 
smoke injected into the channel was increased as the Dean 
number was increased. 

B. CAMERA AND LIGHTING 
1. Radial-Spanwise Plane Photographs 

Radial-spanwise plane photographs were taken using 
generally the same equipment and procedures as described by 
Niver [Ref. 7] and Ligrani and Niver [Ref. 9]. Here, Kodak 
T-Max 400 film was used. The Nikon F3 camera body and 55 
mm, £2.8 lens were used with a Nikon MD-4 motor drive and a 
Nikon MT-2 intervalometer. Best results were obtained with 
an aperature setting of £4.0, a shutter speed of 60, and an 
ASA of 1600. All photographs were taken in sequence with 
the intervalometer set for a pulse time of 0.1 seconds, an 
interval time of 0.1 seconds and in the repeat mode. The 
total time between photographs is 0.2 seconds. A schematic 
of the setup used for photographs of illuminated radial- 


spanwise planes is shown in Figure 27 [Ref. 7]. 
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2. Concave Surface Exterior Photographs 

Photographs were taken from the outside of the 
concave surface at locations encompassing the 50-80 degrees, 
80-110 degrees, and 110-140 degrees streamwise positions. 
The photographs were also taken in sequence using the same 
intervalometer settings as mentioned above. Color Trans 
lighting was used instead of an automatic flash for surface 
illumination. The best results were obtained using Kodak T- 
Max 400 film and ASA set at 1600. The aperature and shutter 
speed were adjusted manually for different lighting 
conditions. Typical settings were f£4.0 and a shutter speed 
ef 30. 

Photographs were taken as the camera viewed one-half 
of the channel width, after adjusting the camera and 
lighting. After a photograph sequence was taken at one Dean 
number, the Dean number was increased and a pause of at 
least sixty seconds was taken in order to allow any 
disturbances from the flow rate adjustment to decay. After 
a range of Dean numbers was photographed at one streamwise 
location, the camera and lighting system were repositioned 
to the next downstream location and the photograph sequence 


was repeated. 
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VIII. RESULTS 


This chapter presents measured mean velocities, mean 
vorticity distributions, power spectra and photographs of flow 
Visualization. All results from this study were obtained for 
Dean numbers from 40 to 230. The three mean velocity 
components and the three components of mean vorticity were 
determined from pressure measurements. All of these pressure 
measurements were taken 120 degrees from the start of 
curvature in a radial-spanwise plane 5.08 cm to 10.16 cm off 
the channel centerline. Spectra of velocity fluctuations were 
obtained at the channel centerline using a subminiature 
crossed hot-wire probe 110 degrees from the start of curvature 
at locations 9.60 and 9.92 cm from the side wall of the 
channel. Flow visualization photographs were obtained at the 
same spanwise-radial plane as the pressure measurements and 
also with the camera viewing the concave surface of the 
channel from the channel exterior. 

This chapter is subdivided into four sections. Part A 
contains a review of previous studies of flow in channels with 
a geometry similar to ours. Part B presents results of five 
hole pressure probe measurements. Part C presents the results 
of hot-wire measurements of velocity fluctuations. Part D 


presents flow visualization photographs taken for this study. 
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A. PREVIOUS STUDIES 

A previous investigation of high aspect ration (40:1) 
curved channel flow with a radius ratio of .979 was made by 
Kelleher et al. [Ref. 5]. These authors measured variations 
of velocity with spanwise direction using a hot-wire at Dean 
numbers from 79 to 113 at a location 135 degrees from the 
start of curvature. Results for different radial positions 
showed significant variations of velocity providing evidence 
of spanwise periodic vortices. 

Finlay et al. [Refs. 1, 8] developed computational 
Simulations of curved channel flow by means of three- 
dimensional integrations of the incompressible Navier-Stokes 
equations. To this author’s knowledge, this is the only 
previously published work which includes a detailed 
description of the structural characteristics of Dean 
vortices. Results from Finlay’s work indicate that 
instabilities in curved channel flow sometimes cause Dean 
vortices to develop with different types of superimposed 
travelling waves. Undulating Dean vortex flow is found 
above a critical Dean number of 35.92, but is not fully 
developed for streamwise distances less than one 
circumference. Due to this, and the fact that travelling 
waves associated with undulating Dean vortex flow have long 
streamwise wavelengths, undulating Dean vortex flow is 
difficult to detect experimentally. As Dean number 


increases, a second type of instability causes twisting Dean 
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vortex flow. Travelling waves associated with twisting Dean 
vortex flow have a much shorter wavelengths than undulating 
Dean vortex flow. Twisting Dean vortices are non-symmetric 
and fully developed for De > 2.5 De,, and streamwise 
distances less than x radians. [Ref. 1] 

Ligrani and Niver [Ref. 9] describe results from a flow 
Visualization study of curved channel flow. These 
experiments were conducted in the same channel as the 
present study. Figure 28 shows a series of photographs 
taken by these authors at the 115 degree streamwise location 
where the concave surface is on the bottom and the convex 
surface is at the top of each photograph. The left hand 
Side of the photographs is 10.16 cm from the channel 
centerline. The width of the area in the photograph spans 
5.08 cm. Flow upwash regions (from concave to convex wall) 
are seen aS narrow "stems" of mushroom shaped smoke 
patterns. The broad "petals" on either side of the upwash 
region contain fluid which is being forced outward and away 
from the convex wall. The first photo shows the initial 
development of a flow instability at De=53. From De=73 to 
De=100 mushroom shaped smoke patterns evidence pairs of 
vortices which increase in radial height with Dean number. 
These vortices are also considered to be symmetric which 
means that the two vortices in each pair are similar. Some 
non-symmetry is evident at De=90 which may be caused by 
partially developed undulating Dean vortex flow or a 
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transient disturbance passing through the channel. At 
De=123 the vortex pairs become non-symmetric. The vortices 
also appear to twist in the spanwise and radial directions. 
This "rocking" motion, characteristic of twisting Dean 
vortex flow, is more noticeable as flow rate increases. 
Eventually, for Dean numbers greater than 192, the mushroom 
shaped patterns evidencing vortex pairs are not 
recognizable. Significant Dean vortex patterns are also 


seen as streamwise location changes. 


B. SURVEYS OF MEAN VELOCITY AND MEAN VORTICITY COMPONENTS 

Each contour plots found in Figures 29 through 58 are 
taken from an 8 by 40 grid of data points spaced .0127 cm 
apart. The Y/d axis in these plots is opposite in direction 
to the R axis previously defined (Figure 1). The Y/d axis 
and the Z/d axis are normalized by the channel height, d, so 
that Y/d = 0 represents the concave surface and Y/d = 1 
represents the convex surface. Z/d 1s measured with respect 
to the channel centerline. 

1. Total Pressure and Streamwise Velocity Distributions 


Total pressure measurements normalized by P are 


mean 


given in Figures 29-31 for De=36-201, where as ley Ue 
average of the maximum and minimum pressure differentials 
measured. A graph of these maximum and minimum pressures 
measured at each Dean number is presented in Figure 32. 


This figure shows that the range of total pressures measured 
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increases nearly linearly with respect to Dean number. 


Streamwise velocity measurements in Figures 33-35 are 


normalized by U, where U is the same mean velocity used to 


calculate the Dean number. 


a. De=49 - De=64 


Figures 29 and 33 show results for De=49 - 64. 


Here, a spatially uniform, undisturbed flow is evident at 


the lower De. As flow rate increases, areas of low pressure 


and low streamwise velocity develop near the concave wall. 


Flow visualization photographs taken at this same location 


(Ref. 9] show a uniform layer of smoke which develops a 


"wavy" appearance along the concave surface as Dean number 


increases. 


b. De=64 - De=86 


Figures 29 and 33 show pressure and velocity 


Variations as the flow rate changes. Here, the areas of low 


pressure and velocity 


height. Gradients of 


nearly extend to the full channel 


pressure and velocity appear to be 


more spanwise periodic at De=73.48 and 86.24 than at 


De=64.49. For Dean numbers in this range, Ligrani and Niver 


indicate the presence 
smoke patterns. These 
the same height. The 
regions correspond to 


concave to the convex 


of spanwise periodic, mushroom shaped 
are symmetric vortices and are not all 
low pressure/low streamwise velocity 

secondary flow upwash regions from the 


wall. 
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c. De=95 - De=137 
For the Dean number range from De=95.52 to 

137.0, Figures 30 and 34 show low pressure, low velocity 
areas which are not spanwise periodic. Flow visualization 
for the same conditions (Figure 28) shows that mushroom 
shaped smoke patterns are full height relative to the 
channel depth. These pairs of vortices are non-symmetric 
and have a slight unsteady "rocking" motion. 

qd. De=137 - De=201 

Figures 31 and 35 show that as Dean number 
increases, time-averaged measurements of pressure and 
velocity indicate that the flow is spanwise uniform. 
Simulations [{Ref. 1] and flow visualization [Ref. 9] 
indicate that vortices are still present for these Dean 
numbers but are much more unsteady, especially in the 
spanwise direction, than at lower Dean numbers. Because of 
the amplitude and period of this unsteadiness, time-averaged 
results show gradients only in the radial direction. 
2. Velocity Perturbation Measurements 

Curved Channel Poiseuille Flow (CCPF) 1s present 
when flow in the curved channel is two-dimensional and 
laminar, without secondary flows. This type of flow exists 
for Dean numbers less than 35.92. CCPF profiles differ from 
straight channel Poiseuille flow velocity profile in that 
the maximum velocity is shifted slightly closer to the 


concave wall. Figures 36 through 38 show plots of the 
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difference between measured streamwise velocity and mean 
velocity from the CCPF solution. The CCPF solution for 
streamwise velocity, U(CCPF), is obtained by determining the 
mass flow rate (and U) from measured data and matching this 
U to the CCPF velocity profile. The data is plotted here, 
using Ug - U(CCPF) normalized by U. This illustrates 
perturbations to CCPF caused by time-averaged secondary 
flows. 

Figure 36 shows that at De=49.27, the measured flow 
is nearly spanwise uniform. As flow rate increases to 
De=53.2 and De=57.5, perturbations in the form of streamwise 
velocity deficits and surpluses are evident. Locations of 
velocity deficits match locations of the low 
pressure/velocity upwash regions shown in Figures 29-35. 

The radial position of the peak velocity deficit at 
De=53.2 (Figure 36) 1s ¥/d=0.3. The Ug > U(GGPR) re=niiee 
from Simulations [Ref. l:pp. 17-18] at De=44.18 show that 
the radial position of the peak deficit is at Y/d=0.25. As 
Dean number increases, our data shows that this radial 
position of the peak velocity deficit moves toward the 
convex wall. At De=137.0 the peak velocity deficit is 
centered at Y/d=-0.5. 

At Dean numbers greater than 137.0 (Figure 38) 
contours of velocity deficits and surplus magnitudes become 


more spanwise uniform. This is consistent with the spanwise 
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uniformity seen in the pressure and velocity data previously 
presented for the same Dean numbers. 

3. Mean Vorticity Distributions 

Streamwise, radial and spanwise components of mean 

vorticity were determined from distributions of the three 
mean velocity components. Figures 39 through 62 present 
results from these vorticity calculations. Vorticity values 
at each point are obtained by numerically differentiating 
velocity values with respect to radial or spanwise distance. 


Streamwise vorticity, wg, radial vorticity, w and spanwise 


igi 


vorticity, w, are given by; 
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respectively. Since measurements were taken in a single 
radial-spanwise plane, velocity changes with respect to 
streamwise distance were not included and assumed to be 
small. In the vorticity contour plots of this section, 
areas of negative vorticity are shown as dashed lines and 
areas of positive vorticity are shown as solid lines. 
a. Streamwise Vorticity 

Streamwise vorticity contour plots are shown in 

Figures 39 and 40. At the lower Dean numbers (53.2-64.49), 


there are no distinct patterns of positive or negative 
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vorticity. Pressure/velocity measurement contour plots 
(Figures 29-35) for Dean numbers between 53.2 and 64.5 show 
disturbances to the flow that appear as spanwise non- 
periodic low pressure/velocity areas near the concave 
surface. As mentioned, flow visualization for these 
conditions shows a wavy smoke layer developing along the 
concave wall. 

Figures 39 and 40 show that, as the flow rate 
increases, areas of positive and negative vorticity are 
present as represented by circular-shaped contour lines. 
From De=95.52 to De=126.3 these vorticity regions are seen 
as concentric circles with higher magnitudes of vorticity 
closer to the center of the circles. At this same Dean 
number range, Figures 30 and 34 show that areas of low 
pressure/velocity extend across the depth of the channel. 
If vorticity and pressure/velocity plots are compared at 
each Dean number, it can be seen that the areas of positive 
and negative vorticity lie on either side of the low 
pressure/velocity regions. 

Figure 40 shows that, as the Dean number 
increases to De=137.0, peak magnitudes of vorticity 
decrease. At De=141.6 and De=146.1, the data shows little 
or no evidence of streamwise vorticity. Here, time-averaged 
pressure and velocity distributions are spanwise uniform due 


to unsteadiness in the flow. 
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To further evaluate the effect of flow rate on 
mean streamwise vorticity measurements, magnitudes of peak 
vorticity are plotted against Dean number. These plots are 
shown in Figures 41-44. Here, a vortex pair is considered 


to be the portion of an wy, contour plot with a negative 


x 
vorticity area just to the left of an area of positive 
vorticity. Typically, three vortex pairs, which will be 
referred to as the left, middle and right pairs, are evident 
in each contour plot. Peak vorticity is then defined as the 
average of the absolute values of peak positive and negative 
vorticity. The left and right vortex pairs (Figures 41 and 
43) show little, if any relationship between peak vorticity 
and Dean number. Figure 42 shows that, for the middle 
vortex pair, peak vorticity values increase as flow rate 
increases for 40 < De < 130. Thus, the strengths of the 
vortices appear to vary across the measurement plane. 

Figure 44 shows the average for all three vortex pairs 
wherein peak vorticity increases slightly with Dean number. 
In all four plots, peak vorticity values decrease above 
De=126.3 due to time-averaging and significant unsteadiness 
in the flow. 


Voreexrcirculation, fF is plotted as a function 


xX! 
of Dean number in Figures 45-48. Streamwise circulation is 
calculated by summing streamwise vorticity values over an 


area Of positive or negative vorticity and multiplying this 


sum by (.0127 cm)“, where .0127 cm is the spacing of the 
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data. The value of lr, plotted 1s the averag-se ae -ae 
absolute values of positive and negative circulation. In 
general, circulation appears to increase with Dean number up 
to De=126.3. Above this Dean number circulation values 
decrease. Figure 48 shows that, when circulation is 
averaged over the entire measured area, it increases with 
Dean number up to De=126.3, after which the time-averaged 
measurements then decrease with Dean number. 
Dawekadlal@vVorkaci ty 

Figures 49 through 51 show contours of radial 
vorticity over the spanwise-radial plane of measurements. 
Since the Ug component of mean velocity is measured with 
better accuracy than other components, radial vorticity 
measurements are more accurate than streamwise vorticity 
measurements. As a result, adjacent positive and negative 
radial vorticity contours provide particularly clear 
indications of Dean vortices. According to Finlay et al. 
(Ref. l:pp. 62-63], the locations of radial vorticity maxima 
and minima might also identify locations of shear layer 
instabilities. 

Figures 49 through 51 show undisturbed flow up 
to De=57.5. At De=57.5 relatively low-level radial 
vorticity magnitudes are present on the right hand side of 
the contour plot. As Dean number increases up to 137.0, 


radial vorticity for each vortex pair increases. At Dean 
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numbers greater than 137.0, no time-averaged radial 
vorticity is evident due to unsteadiness in the flow. 

Figures 52-55 show peak radial vorticity values 
plotted against Dean number. Similar to streamwise peak 
vorticity values, measurements for the left, middle and 
right vortex pairs are plotted separately. Peak radial 
vorticity magnitudes for the entire measured are plotted in 
Figure 55. Here, as before, peak vorticity is defined as 
the average of the absolute values of peak positive and 
negative vorticity. In contrast to peak streamwise 
vorticity data, peak radial vorticity increases with Dean 
number at about the same rate for all three vortex pairs for 
Dean number up to 126. 

c. Spanwise Vorticity 

Figures 56 through 58 show the development of 
spanwise vorticity over the range of Dean numbers where 
measurements were made. Vorticity contours at 
De=36.6 - De=53.2 are spanwise uniform. At higher Dean 
numbers, negative vorticity contours move closer to the 
convex surface at the locations of vortex upwash. Between 
these upwash areas, positive vorticity contours move closer 
to the concave wall. These vorticity contours clearly 
indicate that the location of maximum streamwise velocity is 
close to the convex surface at vortex upwash locations. The 


time-averaged vorticity measurements taken above De=141.6 


5S, 


indicate a spanwise uniform flow. This is consistent with 
the high Dean number measurements previously discussed. 

Figures 59 through 62 show the changes in peak 
Spanwise vorticity magnitudes with Dean number. As with the 
other vorticity components, the three vortex pairs are 
evaluated separately and together. Peak positive and 
negative vorticity magnitudes are plotted separately for 
each vortex pair. Peak negative vorticity values increase 
with Dean number at a faster rate than peak positive 
vorticity. Above De=126.3, peak positive vorticity values 
continue to increase with Dean number while peak negative 
vorticity values decrease. 

4. Detailed Pressure and Vorticity Measurements 
Figures 63 through 74 show detailed plots of 

secondary flow vectors, mean vorticity components, total 
pressure differentials and streamwise velocity. The maximum 
secondary flow vector size is shown in the bottom left hand 
corner of the vector plots. Axes of the plots are 
normalized as previously described. This data was obtained 
by traversing the five hole pressure probe using the same 
methods as for the larger area plots but with thirteen 
radial positions and thirty spanwise positions. The spacing 
between the data points is .076 cm in both radial and 


spanwise directions. 
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a. De=77.9 


Secondary flow vectors at this Dean number 
(Figure 63a) do not show a distinct pattern of secondary 
flow except for the upwash region near the center of the 
Peete VOrticity, Pyo+.,) = Pambient, and streamwise velocity 
data (Figures 63b-65) are similar in quality to previously 
presented plots of data taken at De=73.5 and De=86.2. 

b. De=102.2 

Secondary flow vectors (Figure 66a) show three 
secondary flow upwash regions from concave to convex 
surfaces. Figures 67a,b show that the low pressure/low 
velocity region near the center of the plots coincides with 
an upwash region. Vorticity plots (Figures 68a,b) indicate 
that the vortex pair is symmetric with respect to the upwash 
region. 

c¢. De=119.0 

Secondary flow vectors shown in Figure 69a 
clearly show a vortex pair. Streamwise vorticity contours 
plotted in Figure 69b show that these vortices lie on either 
side of the upwash region. Total pressure and streamwise 
velocity contours (Figure 70) show that the low pressure/ 
velocity region extends from the concave to convex walls. 
Radial and spanwise vorticity plots (Figure 71) indicate 


that this vortex pair is slightly non-symmetric with respect 
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qd. De=131.6 

These plots (Figures 72-74) show upwash regions 
from two vortex pairs. From the secondary flow vectors, it 
to the flow upwash. This data agrees with flow 
visualization taken at De=123 (Figure 28) which shows smoke 
patterns leaning to one side. 
is evident that upwash velocity is much greater than 
downwash velocity. Total pressure, velocity and vorticity 
plots all indicate that the flow disturbance on the right 
Side of the measured area is less intense than the 


disturbance on the left side. 


C. SPECTRA OF HOT-WIRE SIGNAL MEASUREMENTS 

To obtain information on the nature of unsteadiness in 
Dean vortex flow, hot-wire voltage fluctuations resulting 
from velocity fluctuations were analyzed. Spectra obtained 
from Fourier transforms of time-domain signals indicate the 
frequencies of different modes of oscillations in the flow. 
Measurements were taken over the Dean number range of 40 to 
200 using both the subminiature crossed-wire probe and a 
conventional horizontal-wire probe. Measurement techniques 
are described in Chapter 6. 

Oscillations or unsteadiness in the flow contribute to 
energy levels seen in spectral plots. If the flow measured 
is steady, no spectral energy will be apparent. Flow 
oscillations at a specific frequency will result in a 
relatively sharp spectral peak at that frequency and 
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harmonics at multiples of that frequency. If vortices 
oscillate over a range of frequencies, relatively broad 
spectral peaks and their harmonics will be evident. 

The effect of background electronic noise was evaluated 
by collecting data with no flow in the channel. Results 
indicate the presence of sharp spectral peaks with 
principals and harmonics of 60 Hz. This background noise is 
the result of Signals picked up from the power supply and 
Signal processing equipment. These spectral peaks are 
generally four to five orders of magnitude less than 
spectral peak amplitudes with flow. 

In order to select frequency ranges of interest for 
different flow conditions, a frequency range of 0-1000 Hz 
was first examined at all flow conditions. From this 
initial investigation, it was determined that two frequency 
ranges should be used at different flow rates. For Dean 
numbers from 40 to 114, the 0-40 Hz range was used. For 
Dean numbers from 120 to 160, both the 0-40 Hz and 0-300 Hz 
frequency ranges were employed. At Dean numbers from 165 to 
200, the 0-300 Hz range was studied. 

1. De=40 - De=114 

Spectra and time traces for Dean numbers ranging 
from 40 to 114 are presented in Figures 75-89. From De=40 
to De=80, the highest energy level is measured at 2 Hz. 
Additional broad spectral peaks are seen at 8 Hz, 14 Hz, and 


two pairs at 18-22 Hz and 28-32 Hz. Sharp spectral peaks 
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seen in the digital spectral plots at 20 and 40 Hz are 
attributed to digital processing "noise". These same sharp 
peaks are not seen in the real-time spectral plots. Because 
so many events seem to be contributing to spectral energy, 
it is impossible to indicate which peaks are due to 
principal oscillations and which peaks are harmonics. 
Ligrani and Niver [Ref. 9] observed both radial and spanwise 
vortex oscillations in this flow range. The vortices they 
observed generally do not span the full height of the 
channel and vary in height when compared to one another. 

Above De=80, Figures 83-89 show that an additional 
spectral peak appears at about 12 Hz. The amplitude of this 
peak increases as Dean number increases up to De=114. 
Ligrani and Niver have shown that in this same flow range, 
Dean vortices are nearly full channel height and appear to 
be relatively steady. 

2. De=120 - De=160 

Figures 90-101 present spectra and time traces for 
De=120 to De=160 at 0-40 Hz. As flow rate increases, the 
spectral peak that was seen at 12 Hz has decreased in 
amplitude. At De=130, the 28-32 Hz pair of spectral peaks 
has also decreased in amplitude. A single broad peak at 30 
HZ 1S now evident. At De=140 to De=160 there are no 
spectral peaks which stand out relative to background energy 
levels. Subminiature and conventional-sized hot-wires yield 


Similar results for this frequency range. 
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Figures 102-109 present spectra and time traces 
for De=120 to De=160 at 0-300 Hz. Relatively sharp, narrow 
spectral peaks seen at 60, 120, 180 and 240 Hz are due to 
background electronic noise. In this frequency range, 
subminiature and conventional hot-wire sensors yield 
different results. This difference is evident in Figures 
104-105 and Figures 106-107, which show spectra taken at 
De=140 and De=152/153 respectively. Spectra taken when the 
conventional wire is employed show broad spectral peaks not 
evident when the subminiature wire is employed. These 
differences may be due to flow interference or disturbances 
caused by the larger sensor. The differences might also be 
due to the fact that the larger wire senses velocity 
fluctuations over a larger area of the flow. 

Flow visualization observations at this Dean number 
range [{Ref. 9] indicate that Dean vortices are nearly full 
channel height and begin to develop a "rocking" motion. 
Finlay, et al. [Ref. 1, 8] describe this flow phenomena as 
twisting Dean vortex flow. 

3. De=165 - De=200 

Figures 110-120 present spectra and time traces for 
De=165 to De=200. Hot-wire spectra in the De=165 to 200 
range were taken from both subminiature and conventional- 
Sized sensors, and both sensors give similar results. At 
De=165 to De=175, Figures 110-113 show broad principal and 


harmonic spectral peaks at about 60 and 120 Hz. At De=180 
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(Figures 114-115) a second harmonic is evident at about 
180 Hz. These principals and harmonics are believed to 
result from twisting Dean vortex flow. 

Figures 116-120 show that, as the Dean number 
increases above De=184, the principal and harmonics 
previously seen decrease in amplitude. At De=200 no 
spectral peaks are evident. 

Flow visualization observations [{Ref. 9] indicate 
that vortex pairs are twisting and are non-symmetric for 
this Dean number range. Temporal power spectra from 


Simulations of fully developed twisting Dean vortex flow 


[Ref. l:pp. 50-51] show principal and harmonics spaced about 


60 HZ “apare. 


D. FLOW VISUALIZATION 

Flow visualization photographs are presented in Figures 
121 through 134. All photographs for a particular Dean 
number and location were taken in sequence 0.2 seconds 
apart. For the radial-spanwise plane photographs, the top 
of each photo represents the convex wall and the bottom the 
concave wall. The right side is 5.08 cm off the channel 
centerline. For the exterior view photographs, the left 
Side of the photos represent the channel centerline and the 
right side is the channel side-wall. Streamwise distance 
increases from top to bottom of each exterior view 


photograph. 
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1. Radial-Spanwise Plane 

Figure 121 presents four photograph sequences taken 
in the radial-spanwise plane. The white smoke was injected 
close to the bottom of the channel inlet and remains along 
the bottom of the 2.44 m straight section as it enters the 
test section. In the curved section, the smoke remains 
along the convex surface until it is disturbed. Because of 
this method of injecting smoke, flow disturbances 
originating along the convex surface contain much denser 
smoke than disturbances originating near the concave surface 
Meer. 9:p. 13). 

Figure 12la presents a sequence of five photographs 
taken at De=51.9. These photographs show a disturbance of 
the flow along the concave surface. This wavy smoke layer 
is neither spanwise periodic or steady. The disturbances 
appear to move in both the spanwise and radial directions. 

Figure 121b presents a sequence of four photographs 
taken at De=64.4. These photographs show mushroom-shaped 
smoke patterns which oscillate in the radial and spanwise 
directions. The smoke patterns are, in addition, not all 
the same height. 

Figure 121c presents a sequence of four photographs 
taken at De=95.4. Here, the mushroom-shaped patterns are 


Spanwise periodic and relatively stable. When individual 
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smoke patterns are distorted, they all appear to lean or 
distort in the same spanwise direction. 

Figure 121d presents a sequence of three photographs 
taken at De=107.7. Here, the mushroom-shaped smoke patterns 
indicate vortices that are nearly full channel height. 
Distortion in the smoke patterns indicates unsteadiness in 
the flow, possibly due to a transient disturbance. 

2. Exterior View - Concave Surface 

Figures 122 through 134 present photographs of a 
portion of the exterior concave surface of the channel. 
These photographs are shown in the sequence they were taken, 
clockwise from upper left. The secondary flow evident in 
these photographs appears as tubes of smoke. These tubes of 
smoke surround vortex pairs which are relatively void of 
smoke (Ref. 9:pp. 22-23}. These photographs do not show 
clear evidence of any type of travelling wave superimposed 
on the flow. It 1s evident however, that the tubes of smoke 
are not always Spanwise periodic. In many sequences, an 
apparent merging or splitting of vortex pairs is seen, 
probably as a result of vortex pair destruction and 
creation. Shortly after this activity, the vortex pairs 


tend to retain their "normal" spanwise spacing. 
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IX. SUMMARY AND CONCLUSIONS 


The development and structural characteristics of Dean 
vortices in a curved channel were measured and evaluated. 
Methods employed for experiments include the use of a 
miniature five hole pressure probe, a subminiature crossed 
hot-wire probe and flow visualization photography. 

The five hole pressure probe was developed specifically 
for use in the curved channel. Its qualification included 
calibration in a wind tunnel and the development of spatial 
resolution and induced downwash velocity correction schemes. 
Pressure measurements were used to determine the three mean 
velocity components of Dean vortex flow. From these, the 
three components of mean vorticity were also determined. 

Pressure measurements show that, at low Dean numbers, 
spanwise non-periodic disturbances develop near the concave 
surface of the channel. Above De=57.5, vortices were 
measured and vorticity levels calculated. It was determined 
that vortex pairs tend to become spanwise non-periodic as 
Dean numbers increases above De=95. Levels of vorticity 
increase as Dean number increases, until steady vortices 
caused by twisting Dean vortex flow are present. Because 
the length scales of the spanwise unsteadiness are greater 
than spanwise vortex spacing for De > 140, pressure 


measurements show gradients only in the radial direction. 
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Hence, at these Dean numbers pressure measurements do not 
resolve local spatial variations due to twisting Dean vortex 
flow. 

Velocity fluctuations caused by unsteadiness in Dean 
vortex flow were measured with a subminiature crossed-wire 
hot-wire as well as with a conventional horizontal hot-wire. 
Spectral analysis of these fluctuations indicate that, at 
relatively low Dean numbers (De=40 - De=114), several low 
frequency modes of vortex oscillation are present. For 
De < 90-100 these are associated with partially developed 
Dean vortex flow. As Dean number increases greater than 
160, spectra show a principal peak and harmonic peaks, 
evidencing the presence of twisting Dean vortices. These 
measurements above De=160 compare favorably with the results 
of simulations conducted using similar flow conditions. 

Dean vortex flow was observed using flow visualization 
techniques. The characteristics of pairs of Dean vortices 
appears to be strongly affected by the local small 
disturbances which initiate vortex development. When smoke 
patterns are viewed from the outside of the concave surface, 
they appear as tubes of smoke, each of which surrounds a 
vortex pair. Because of new development and destruction of 
vortex pairs, these occasionally appear to merge and split. 
Before and after such activity, a specific spanwise spacing 


seems to be maintained. 
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APPENDIX A 
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Figure 1. Curved Channel Geometry 


Figure 2. Schematic of One Dean Vortex Pair 
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Figure 9. Pressure Probe Calibration Sled 
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Figure 11. Cp Pitch vs. Pitch, 1.9 m/s 
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Figure 12. Ch Total vs. Pitch, 1.9 m/s 





Figure 13. Cp Total-Static vs. Pitch, 1.9 m/s 
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Figure 19. Successive Positions of the Five Hole Pressure 
Probe and Probe Ports, Looking Downstream 
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Figure 24. Secondary Velocity Vectors, De=100 
a) Spatial Resolution Correction Only 
b) Spatial Resolution and Downwash Velocity 
Correction 
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Figure 26. a) Hot-wire Measurement Setup 


b) Oscilliscope, Bridge and Filter 
Amplifier 
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Figure 27. Radial-Spanwise Plane Flow Visualization 
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Figure 28. Radial-Spanwise Plane Flow Visualization 
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Figure 29. Total Pressure Contours, De=49.3 to De=95.5 
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Figure 31. Total Pressure Contours, De=126.3 to De=201.8 
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Figure 32. Total Pressure Differentials 
Maximum and Minimum 
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Figure 33. Streamwise Velocity Contours, 
De=49.3 to De=95.5 
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Figure 35. Streamwise Velocity Contours, 
De=126.3 to De=201.8 
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Figure 36. Velocity Perturbation Contours, 
De=49.3 to De=95.5 
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Figure 38. Velocity Perturbation Contours, 
De=126.3 to De=201.8 
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Figure 39. Streamwise Vorticity Contours, 
De=53.2 to De=100.3 
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Figure 40. Streamwise Vorticity Contours, 
De=100.3 to De=146.1 
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Figure 45. Vortex Circulation vs. Dean #, 
Left Vortex Pair 
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Figure 46. Vortex Circulation vs. Dean #, 
Middle Vortex Pair 
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Figure 47. Vortex Circulation vs. Dean #, 
Right Vortex Pair 
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Figure 48. Vortex Circulation vs. Dean #, 
Average Vortex Pair 
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Figure 49. Radial Vorticity Contours, 
De=36.6 to De=86.2 
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Figure 51. Radial Vorticity Contours, 
De=137.0 to De=230.7 
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Figure 54. Peak Radial Vorticity vs. Dean #, 
Right Vortex Pair 
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Figure 56. Spanwise Vorticity Contours, 
De=36.6 to De=86.2 
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Figure 57. Spanwise Vorticity Contours, 
De=86.2 to De=137.0 


izes 


DERN #3236524 

















































OMEaGRz Gl 7S) a 8 28 3 8 a 666 68 
RANGES 
@:2-338 <-282 
1 :>-2@6 <-138 
2 :>-138 <-82 
3 :>-8B <-42 
4 :9-48 < 48 
S :> 40 < 82 
6 :> @@ < 138 
? :> 138 < 228 
8 :> 208 < 272 
DEAN @# 176.1 
? ? ? ? ? ? ? as 8 
ee a ae oe el ee ee er 
DERN @# 155.3 
S 6 a 7 oe 7 ae ? 777 ? 
<=> yp Trt yT=s =pS=>= Se ae 
DEAN @ 146.1 
72 ae 7 ? ? ? 
DERN # 141.6 
? ? S ? ? ? ? ? ? ? 
eee eeu 
g==ET=aLSr TRESS Es=e==eT=FS- 2. _-3~RI=E: 
DEAN # 137 
.8 ? ? ? aa 
~, — —— 
se) ee = ee 
4 = VG SE —~——_ Pew a oe 
e Sg S Sse FSS FEST SLSss=22ZE>S-~--27- ==: 
",) 


Figure 58. Spanwise Vorticity Contours, 
De=137.0 to De=230.7 
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Figure 63. De=77.9, a) Secondary Flow Vectors 
b) Streamwise Vorticity Contours 
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Figure 65. De=77.9, a) Radial Vorticity Contours 
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Figure 66. De=102.2, a) Secondary Flow Vectors 
b) Streamwise Vorticity Contours 
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Figure 67. De=102.2, a) Total Pressure Contours 
b) Streamwise Velocity Contours 
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Figure 88. De=110, a) Subminiature Hot-wire Spectra, 
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b) Sample Time Trace 
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Figure 94. De=130, a) Subminiature Hot-wire Spectra 
Digital Analysis 
b) Sample Time Trace 


160 


40. 000 


Serta 


AF .0500 
(HZ) 


» 
Pipe UGE ey 


oer 
=. 382 
fei tt | i. 


Tt MN Ce CD 
as: a | 
< al Lil 
“ <= 
=> 


Figure 95. De=134, Subminiature Hot-wire Spectra, 
Real Time Analysis 


161 


40. 000 


SJ 
xe 
Ce 
i) 
nn -_— 
O wn 
- ae 
Li 0lti<‘ 
<j 
> 
eS) 
PLL 
Lu) 
> 
Ss 
2 
< 
or 
Se 
Pa 
= 
=D 
ere 
| 
rs 





CA 
CCT. 


uJ 
<t 


V? /HZ 
LOG 
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Figure 105. De=140, a) Conventional Hot-wire Spectra, 
Digital Analysis 
b) Sample Time Trace 
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Figure 109. De=160, a) Subminiature Hot-wire Spectra, 
Digital Analysis 
b) Sample Time Trace 
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Digital Analysis 
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Figure 112. De=170, Subminiature Hot-wire Spectra, 
Real Time Analysis 
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Figure 113. De=175, a) Conventional Hot-wire Spectra, 
Digital Analysis 
b) Sample Time Trace 
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De=180, a) Conventional Hot-wire Spectra, 
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Figure 116. De=184, Subminiature Hot-wire Spectra, 
Real Time Analysis 
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Figure 117. De=190, a) Subminiature Hot-wire Spectra, 
Digital Analysis 
b) Sample Time Trace 
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Figure 118. De=190, a) Conventional Hot-wire Spectra, 
Digital Analysis 
b) Sample Time Trace 
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Figure 119. De=190, Subminiature Hot-wire Spectra, 
Real Time Analysis 
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Visualization 
a) De=51.9, b) De=64.4, c) De=95.4 da) De=107.7 


187 


SS 


. 


SS 


Sektenatae” 
nnn > 
SHOR Hee 


SS 


ANS 


SS 


AQAA’ 


oN Ny *, ‘aaa ta es 
SENNA Seta ate 
. s 


es 


SS 





ie 


LireZa 


e 


1sua 


122. Concave Surface Flow V 
De=88.5, 


igure 


EK 


Lone 


=o Um 


50° 


188 


————E— 


Figure 123. 


SS 


< SAAN 


Concave Surface Flow 
De=97.0, 50° -80° 


189 





Visualization, 


ASI 


Soy 


Ny 


SQ 


as 


WHY RBA GG OS 


TSN 





Figure 124. Concave Surface Flow Visualization, 
De=106.6, 50° -80° 
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Figure 127. Concave Surface Flow Visualization, 
De=165.6, 50° -80° 
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Figure 129. Concave Surface Flow Visualization, 


De=64.5, 80°-110° 
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Figure 130. Concave Surface Flow Visualization, 
De=74.8, 80°-110° 
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Figure 132. Concave Surface Flow 


Visualization, 
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Figure 135. Uncertainty Analysis Flow Chart 
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APPENDIX B 


UNCERTAINTY ANALYSIS 
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I. DEAN NUMBER UNCERTAINTY 


Uncertainty in Dean number calculation is derived from 
the uncertainty in the measurement of the pressure drop 
across an ASME standard orifice. A sample set of 


uncertainty quantities for De=150 is listed below; 


A. CELESCO TRANSDUCER CALIBRATION UNCERTAINTY 


OEP se KCC): sale H50 


manometer ~ — 


§E = + .0007 Volts 


oC 


#20012 in. H,0/Volt 


B. PRESSURE DROP UNCERTAINTY 
$E = + .047 Volts 


SAP =+ .005 in. H50 


C. MASS FLOW RATE UNCERTAINTY 
6A. = + .012 in.? 
§K = + .001 
§Y = + .02 


.002 lbm/ft? 


ray 
o 
Il 

it 


mH 

| 
tl 

Tae 


.0004 lbm/sec (95% confidence level) 
§Ach = + .0052 in.? 


6De = + 4.4 (95% confidence level) 
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II. FIVE HOLE PROBE VELOCITY MEASUREMENT UNCERTAINTY 


Uncertainty in five hole probe velocity components is 
derived from the uncertainty in actual pressure measurements 
as well as the inherent uncertainty in the probe 
calibration. The flow chart used in determining the 
uncertainty of velocity components is shown in Figure 135. 

A sample set of uncertainty quantities for De=100.3, Y/d=.5, 
Z/a=6.0 is listed below. Numbers refer to the position in 
the flow chart. 

A. TRANSDUCER CALIBRATION COEFFICIENT UNCERTAINTY, 

VALYDINE TRANSDUCERS 


6 P = Ft 01 H50 


manometer 


§E = + .0007 Volts 


te 


§c = -.000072 in. H20/Volt (bias error is disregarded 
Since it will cancel out when pressure differentials are 


calculated). 
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B. PROBE CALIBRATION UNCERTAINTY, PITCH = 0°, YAW = -10° 
§E = + .0004 Volts 
6 Py =. OOO6 lin. H50 
6 Po te OOO lim. H50 
1) 6 P23 = te EO. 5 sah. H50 
6 Py = + .0002 in. H50 
SIs = + .0003 1Nn. H.0O 
je 2— see OOO 04am . H50 


§P, = + .0007 in. H,0 


§P = + .00015 in. H50 

6 C = + .025 

2) py 
= + 

6 Con + . OOF) 

| §Cpts = + -045 
6K, = + .004 degrees" 1 
S)) 6 Ky, = + .004 degrees 1 
al 


Gee See GmCegimees © 
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C. CURVED CHANNEL MEASUREMENT UNCERTAINTY 


SE. §C as before 


6 Po =) tt OO Oza Line H50 
4) §P, = + .0002 in. H,0 
§P, = + .0002 in. H30 
6P; = + .0002 in. H50 
§P = + .0001 in. H,0 
SE @ = + .134 
5) py 
= 
6 Cop te eee tals 
§Yaw = + 1.33 degrees 
6) 


§ Pitch = + 1.4 degrees 


7) 6 (Protal = Potatic) = + .0004 in. H50 


8) 6U = + .063 m/s 


6U,g, = + .062 m/s 
=) §6U, = + .062 m/s 
6U, = + .030 m/s 
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APPENDIX C 


SOFTWARE DIRECTORY 


This Appendix gives a listing of the various programs 
used in the thesis. Each program listing contains a summary 
of how the program is used, user inputs, program outputs and 
additional features, if any. Programs are listed by type; 
five hole, pressure probe, hot-wire spectra or plotting 
routines. All programs were written in BASIC 4.0 for use on 


the HP9836S computer. 
I. FIVE HOLE PRESSURE PROBE PROGRAMS 


A. ORIENT: 
Program to orient five hole probe in the calibration 
wind tunnel before calibration begins. 


user inputs: coefficients of transducers calibration 


program output: Pi, P2, P3, P4, P5, P, Coy Sop 


additional features: allows user to continue through the 
program and adjust the probe calibration sled until 


satisfactory Coy and Cop are attained 


On, GC = ate 


(Coy = pp 
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B. PROCAL: 


Program used for calibration of the five hole pressure 


probe in the calibration wind tunnel. 


user inputs: a) 
b) 
Cc) 


d) 


program outputs: 


number of yaw points. 
number of pitch points. 
pitch and yaw resolution. 
pressure offset corrections due to 
pitching of transducers. 
a) data file of calibration coefficients; 
pitch, yaw, Coy’ Cop: Cote Cos: Cots: 


b) data file of pressures; pitch, yaw, 


Par Por P3r Par P5r Pgtaticr Ptotal: 


additional features: a) allows user to calibrate 


transducers or input calibration 
coefficients. 

b) allows to vary yaw at constant 
pitch or vary pitch at constant 
yaw. 

c) can accommodate five, six or seven 
transducers. 

ad) allows for simultaneous printout 


of data if desired. 
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C. ALIGN: 

Program to align the five hole pressure probe in the 
curved channel at zero pitch and zero yaw angles. 
Interfaces with MITAS controller to position probe at a 
desired location, repositions probe to original position 
upon completion of the program. 
user input: a) original position of the probe. 

b) desired position of the probe. 


program output: total velocity, Py, Po, P3, 


“pp: “py: 
ra Pc. 
additional features: allows user to re-align the probe 


between measurement sequences. 


D. FIVE HOLE: 

Program for control of five hole probe positioning and 
data acquisition in the curved channel. Interfaces with 
MITAS controller to control the probe position. 
user inputs: a) initial position of pressure probe. 

b) number of x and y points to take data at. 

c) resolution of data. 

ad) ambient pressure. 

e) delay time to start the experiment. 
program outputs: a) data file of raw pressures; 


Set eo Ean ay Ps, P, Dean no. 


Z09 


b) data file of uncorrected pressure 
coefficients; X, Y, Choy: Cop 
additional features: a) reads temperature and calculates 
ambient density and Dean number for 
each group of pressure measurements. 
b) averages Dean number at end of 
experiment. 
c) allows user to calibrate trans- 


ducers or input calibration 


coefficients. 


E. PADJUST: 
Program to implement spatial resolution correction as 
described in Chapter IV. 
user inputs: a) distance between port 1 and center of other 
pressure ports (YLEN, ZLEN), .016 in. for 
miniature five hole pressure probe. 
b) data file and name for raw pressure data, 
same raw pressure data as is output from 
Riv i HOLE. 


program output: corrected file of pressure coefficients. 


F. VELOCITY: 

Program to calculate pitch and yaw angles for curved 
channel flow meaSurements. From this data the three mean 
velocity components at each measurement location are 
determined. 
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user inputs: a) data file of calibration coefficients and 
pressure data as output from program 
PROCAL. Program reads in this data for 
linear/bilinear interpolation. 
b) values of A (DELY, DELZ) for y and z 
direction downwash velocity correction 
(Chapter IV). Typically a value of 0.1 
was used for A. 
c) data file of pressure coefficients as output 
by program PADJUST. 
program output: data file of pressures and velocities for 


each measured location as follows; 


Y, 2, P total Us Ugr Uy, Uz: 


G. CCPF: 

Program to calculate Curved Channel Poiseuille Flow at a 
given Dean number. The value of CCPF at each radial 
location is subtracted from actual streamwise velocity 
values at each measured location. 
user input: data file of velocities as output by program 

velocity. 
program output: U, - UCCPF values as follows; 


Y,Z,UCCPF,U,,-UCCPF 


Zee 


H. VORTICITY: 

Program to calculate the three mean components of 
vorticity at each measured location. Methods used for 
calculating vorticity are presented in Chapter VIII. 
user input: data file of velocities as output by program 

VELOCITY. 
program output: data file vorticities as follows; 


Y Z A) WO + 03 


Ie CIRC: 

Program to calculate peak vorticity levels and 
circulation values from a particular set of vorticity data. 
user input: a) data file of vorticities as output by program 

VORTPCiDY. 

b) threshold limits, as a percentage, for 
omitting data points when circulation is 
calculated. For example, a 20% threshold 
will cause the program to omit all vorticity 
values less than 20% of peak vorticity. 

c) limits of Y and Z to "look at" when 
calculating circulation or peak vorticities. 
This allows the user to calculate circulation 
for a particular pair of vortices as "seen!" 


from previous data plots. 
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program output: hard copy and/or binary data file of Y and Z 
limits, peak vorticities and circulation for 


a particular vortex pair. 


Program to sample orifice pressure drop and calculate 
Dean number (see Chapter II). 
user input: a) ambient pressure 
b) transducer calibration coefficient, or allows 
user to calibrate transducer. 

program output: a) mass flow rate 

b) orifice pressure drop 

c) Dean number 
additional features: allows user to adjust throttle valve 

(flow rate) and then calculates 


another Dean number. 


K. GET-DATA: 
Program to access various data files and printout 
applicable data. 
user input: a) type of data file (Pressures, Pressure 
coefficients, Velocity). 
b) Y and Z location to sample. 


program output: list of data for sample point. 
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II. HOT=-WIRE VOLTAGE SIGNAL, DATA ACQUISITION AND 
PROCESSING 
A. HOT-WIRE: 

Program to interface with HP67755A high speed data 
acquisition system and collect a sample set of hot-wire 
voltage signals. 
user input: a) Dean number 

b) Sampling frequency 


program output: data file of voltage values. 


B. CONTOUR: 

General contour plotting routine has been adapted for 
specific plotting purposes. Plots contour lines between 
predetermined levels. 
user inputs: a) type of data file to plot (pressure, 

velocity, vorticity, ere. 


b) bytes of data to skip over. 


c) number of data files to average 
(maximum is 20). 

ad) option to plot average values, 
average of deviations from mean value, 
or average of standard deviations from 


mean value. 
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e) number of contour levels to plot. 
f) allows for user to define levels or will 
equispace levels. 
program output: contour plot and table of levels, labels 


local extrema. 


C. QPLOT: 
General X-Y plotting routine. 
user inputs: a) number of columns in data to be read. 
b) column of independent variable (x). 
c) column of dependent variable (y). 
d) arrangement of data file. 
e) titles, subtitles, etc. 
program output: an x-y plot of data as input can include 


up to five curves per plot. 


D. O-SCOPE: 
Adaptation of QPLOT to plot hot-wire voltage signal vs. 
time. 
user inputs: a) time to begin plotting (point number) 
b) total time range to plot (determined from 
frequency of data acquisition). 


program output: hot-wire voltage vs. time plot. 
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E. SAMPLTCHAN: 

Program to sample data from up to seven channels from 
HP3497A data acquisition system. Used for plotting pressure 
signal vs. time. 
user inputs: a) number of channels to sample. 

b) number of samples to take per channel. 
c) time between samples. 
program output: signal vs. time plots for up to seven 


channels. 


F. VECTOR: 
Program to plot secondary flow vectors ina flow field, 
largest vector is labeled (5% of subplot width). 
user inputs: velocity data file (output from program 
VELOCITY. 


program outputs: secondary velocity vectors. 


G. PLTRAV: 

Program to plot traverses of pressure probe measurements 
will plot up to ten radial traverses, for different spanwise 
locations ina "three-dimensional" format. 
user inputs: a) velocity data file (output form program 

VELOCITY) . 
b) Spanwise positions to plot. 


program outputs: radial traverse plots. 
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H. 3DPLOT: 


Program to plot seven radial-spanwise contour plots from 
curved channel measurements (See Figure 33). 
user inputs: a) type of plot desired, vorticity, total 
pressure, streamwise velocity or U,-UCCPF. 
b) Appropriate data files. 


program output: Contour plots. 
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